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1. Introduction 

The notion of entanglement was introduced by Schrodinger 1 1 1 who considered as the 
essential feature of quantum mechanics the fact that when two particles interact by a known 
force and then separate, they can no longer be considered as independent. This ought to be 
evident in the case of a hydrogenic atom for which, by some mechanism, the interparticle 
interaction is lost and the atom is ionized. The entangled nature of the state must be evident in 
the structure of the two-particle wave function, but this does not appear to have been examined 
in detail. The aim of this note is to study this situation for an exactly solvable two-particle 
system. 

A suitable one dimensional model was introduced in 1 2 1 to describe the interaction of 
a hydrogenic atom with a "metal" surface, such that once the atom penetrates the metal the 
nuclear "charge" is screened to zero. It was found that the problem could be re-expressed as 
the Wiener-Hopf problem introduced in 1947 to describe the reflection of an electromagnetic 
wave from a linear coastline solved in 1952 by Bazer and Karp. Their work is reviewed in 
and from the results, a formal expression for the two-body wave function was obtained along 
with an exact formula for the reflection coefficient as a function of the incident energy E. 

In the atomic center of mass system (total mass M, center of mass coordinate R, reduced 
mass /i and relative coordinate r) the interaction potential for the model is V(R, r) = —X5(r) 
outside the metal (R > 0) and inside {R < 0). In terms of reduced variables: 



the wave function is (the sign in the exponent of the Fourier transform has been altered from 
that in 1 2 1 to reflect standard usage) 





ip(R,y) = exp[-iKR-a\y\} + 



a 2 [a+{K)} 2 



cxp[iKR — a\y\] 



2K 2 
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aa+(K) / cxphifcfl-MV^g] „ +(jfe) dfc (fl > Q)> (L2) 



27ri J (fc + X) [a - Vfc 2 - fcg] 

m^n^-ivi^gi (L3) 



2tti / (fc + JQ-y/fc 2 - fcg 
where, cr ± (fc) are the Wiener-Hopf factors |4| of the kernel 

a(fc) = 1 _^^ = ^) (1 . 4 ) 

and the contours for the integrals in II. 2\ and II. 3> surround branch cuts associated with ±fco 
(which, for technical reasons, have infinitesimal imaginary parts) in the upper and lower half 
planes. 

The leading term on the right hand side of II. 2> represents an incoming atom and is 
normalized to unit amplitude; the second term represents the atom reflected from the surface 
with reflection coefficient 

» = ^. 

The task here is to analyze the integrals in II. 21 and II. 3> which describe the entanglement 
of the atomic particles in the "ionized" state. Since for R < there is no direct interaction 
between the two particles comprising the atom, a dependence of the wave function on the 
relative coordinate y indicates the quantum entanglement of the two particles. 

In the following section we change II. 3> and ll.4> into somewhat more convenient forms 
and express all the components of the Wiener-Hopf solution in terms of standard functions. 
In the concluding section details of the wave functions are presented and we find that indeed, 
even for R large and negative, the two-particle wave function is not separable. 

2. Calculation 

We begin by investigating the Wiener-Hopf factorization of er(fc) II. 41 . Only the magnitude 
of the + factor at the value k = K was needed in 1 2 1 ; we now require its complex value for 
Im k > 0. The calculation is simplified somewhat by introducing the modified kernel 

in terms of which II. 2\ and (II .31 become 

1>(R, y) = exp[-iKR - a\y\] + 2 (|— ^ S 2 + (K) exp[iKR - a\y\] - *(R, y) (2.2) 



® {R ' y) ~ YTk- Q S+i ] T m Vfc^ — s + {k){K*-k*) (i?>0) ' (23) 

where the contour encloses the branch cut running from — ko in the lower half plane, parallel 
to the real axis to the imaginary axis and then to — ioo. For R < 0, 



2aKS+(K) I dfc [kTh exp[-ikR~\y\^k^kj} 
Wy) = K + ko f 2^i yj—Vo S + (k)(K*-k>) (2 ' 4) 

where the contour surrounds the branch cut running from ko, parallel to the real axis up to 
the imaginary axis and then to iac. An alternative derivation of d2.2i — d2~4l > is outlined in 
Appendix A. 
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Now Elm, 



S+(k) = exp[-J(fc)] 



1 f c 



Log 



1 



fort 



du 



u - 



(2.5) 



where the contour runs along the real axis indented below u = fc. By breaking the range into 
[— oo, 0] U [0, oo] and combining the two integrals we have 



J(fc) 



Log 



1 



du 



u 2 — fc 2 



(2.6) 



Next, by letting u — > — is and, since there are now no singularities in the first quadrant, 
rotating the contour back to the real axis, with s — ► fcs, we find 



J(fc) = - 

7T 



S 2 + 1 



Log 



1 



V* 2 + (Vfc) 2 



ds 



-LogWs 2 + (fco/fc) 2 - ia/k] 



1_ Log[s 2 + (fcp/fc) 2 
2tt y s 2 + 1 



ds. (2.7) 



The second integral, as can be found in tables, is 7rLog[l + (fco/fc)]; the first integral is 
evaluated in Appendix B yielding 



S+(k) 



fc + fco 



K 



exp 



2i 



Im Tio 



y/fc 2 - fc 2 + k 
K + k 



(2.8) 



where T'\2{z) is the Arctangent integral function, defined in JB.3i . Formula (12. 81 is valid for 
Im fc > 0, except for fc = ±K, where there is a confluence of singularities. However, this 
case can be extracted from the results given in 1 2 1 whence we find 



+ ( ' = \^2K~ 6XP 
S+(K)S + (-K) = i. 



— [U 2 (-a/K)-U 2 (a/K) 
Zir 



(2.9) 
(2.10) 



3. Results and conclusions 

We turn now to the wave function &2.4i in the interaction-free region. At this point the 
imaginary part of fco can be set to zero and on the part of the contour surrounding the 
vertical portion of the branch cut, where the factor exp(— ifci?) is a decaying exponential, 
the contribution to the integral will be small and we shall neglect it. Similarly, since the 
contribution of the contour in the upper half plane about the interval [0, fco] is exponentially 
small compared to the part in the lower half plane, for R < we have 



m y)**«S + (K) /^^g^-gv^l, (3 . 1} 



K + k a +v 'J 27rVfco-fc S+(k)(K 2 - k 2 ) 
Note that, from i\.5\ . as fco — ► 0, 1Z = 1 signifying total reflection. Accordingly, from 
( 13. It . \ip(R, y)\ 2 =0 apart from a rapidly decaying evanescent wave due to the neglected part 
of the contour. 

We have evaluated ( I3.lt numerically for a = 1, fco = % ar, d present the square of its 
amplitude as a function of |y| and \R\ for a = 1 and fco = 2 in Fig. 1. In Fig. 2 we show the 
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Figure 1. Square of the wave function amplitude \i))(R, y)\ 2 ' a = 1, ko = 2 



real and imaginary parts of ip(R, y). For R large and negative, but \y\ <C \R\, the integral can 
be estimated by setting k = in the integrand, apart from the factor exp(— ikR), leading to 

MR, y) » ; a e -i[*olvl+*-] (3.2) 

mK*Ry/2ko(K + k ) 

Im Tio 



7T V K 

Thus, the amplitude decays as \R\~ 1 and even far away from the point where the particles 
are decoupled, the relative coordinate is present in the phase. For larger values of \y\, the 
amplitude appears to be o(|y| -3 / 4 ) and oscillates as y — > oo, which is illustrated in Fig. 3 
showing — 10, y)\ 2 for a = 1, ko = 2. 

A natural quantity to examine is the expected value 

_ J^lyMR^dy 



of the relative particle displacement . Our numerical integration of the numerator in( 13. 3> 
appears to diverge in the interaction-free region R < 0. To obtain a reliable result will require 
a more detailed asymptotic study of the wave function for \y\/R — * oo. Finally, in Figure 4 
we show the square amplitude of the wave function il.2\ for R > with a — 1 and fco = 2 
(including the incident wave). 

For R S> 1 the integral in dl .31 has been treated by Kay 1 3 1 who used the method of 
steepest descent. From his results we find for the integral in (12. 3t . 



1 ' yJ 47r 1 /2(fc oi? )l/2 a 2 + K 2^2 ^ KK+ 1 ^2 V") 

1 U 2 (a 2 /K 2 )~2U 2 (a/K)-lmTi 1 ^ + ia 

= y/R < 1 . This ref 
representing the atomic state. 



" + 2 zy 1 ' zw ' \ K -k + ±k e 

where f = y/R < 1. This represents an outgoing wave with no trace of the factor exp(—a|y|) 
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0.2 




Figure 2. Real (left) and Imaginary (right) parts of ip{R, y): a = 1, kg = 2 



l^(-10,2/)| ; 
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Figure 3. 10, y)| 2 : a = 1, ko = 2 
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where L = d R + d 2 + k 2 . and has the Green function G(R — R',y — y';ko), which in 
momentum space is simply [P^ + p 2 — k 2 ] -1 



oo poo 



Figure 4. \ip(R, 0)| 2 (upper curve) and \ip(R, .5)| 2 for R > 0: a = 1, k = 2 

Appendix A. 

The two-particle Schrodinger equation is 

Lip(R,y) = V(R,y)ip(R,y) (A.l) 

r + d 2 y + k% 
ice is simply [j 
The solution of ( IA.H is 

fl>{R,y)= I I G(R - R'y - y';k )V(R',y')ip(R',y') dR' dy'. 

J — oo J — OO 

In terms of h(R) = i^(R, 0) this reduces to 

p OO 

rp(R,y)=2a I G(R — R', y; k )h(R') dR' (A.2) 
Jo 

and, in particular, 

/oo 
G(R- R',0;k a )h + (R')dR' (A.3) 
-OO 

where h+(R) = h(R)Q{R), h-(R) = h(R)G(-R). The Fourier transforms of these 
functions H + (k) and H_(k) are analytic in the upper and lower half fc-plane, respectively. 
Since the integral in iA.'H is a convolution, by taking the Fourier transform ( IA.3> becomes 

H+{k)+H-{k) = 2ag a (k)H + (k) (A.4) 

where go(k) = \j2\Jk? — k 2 is the Fourier transform of G(R, 0; fco). Therefore, 

H-{k) k 2 -K 2 



H+{k) k 2 -k, 



2 ^2 S(k) (A.5) 



n 



where 
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and S+(k)(S-(k)) is analytic and free of zeros for Im fc > (Im fc < 0). Equation iA.bt 
implies that (fc 2 — K 2 )S + (k)H + (fc) /(fc + fc ) is a bounded entire function and is therefore a 
constant a. Consequently, 

/>oo+ie i | / 

= g *t* ^ e-^dfc (A.7) 



y_ 00+ic s'+(fc)(fc a -Jf 2 ) 

with Im _K" < c < Im fcp. Thus, we have from iA.H , for — oo < R < oo, 



\ aa A 004 "* „ /fco + fc e - ifc «-|ylV fe2 - fc o 
^ = W^** V A^fc (fc 2 -^)5 + (fc) - (A ' 8) 
For R > the contour can be closed into the lower half plane, which contains the two poles 
±K and the branch point — fco yielding i2.2i - $23l . while for R < the contour can be closed 
into the upper half plane yielding &2A\ . 



Appendix B. 

Here we evaluate the integral 



Log[Vx 2 + c 2 + d] 

5— dx (B.l) 

needed to obtain S+(fc). We begin indirectly by setting 

^ = ycot(e/2) B = ytan(fl/2) 
1 + vV + 1 1 + Vl + 2/ 2 ' 



Then 

and 
But, 



Therefore, 



where 



A + B ra\ 



du 



1 - AB 

tan _1 (ycsc(6))) = tan -1 (A) +tan _1 (S). 

csc(0)tan _1 (ucsc(0)) = / — y ' 

//.n 1, d f A tan _1 (u) , 

csc(0) tan" 1 (A) = — - / ^ du 

do J u 

1 d f B tan _1 (w) , 

csc(0)tan _1 (i3) = — / ^ du. 

d9 J u 

p>?^M d „ =Ti3(A) _ Ti2(B)+c , 

Ji Vti -1 



, tan (u) , „ 
Ti 2 (z)= / ^du (B.3) 



is the Arctangent integral function and C is independent of 0. However, for 9 — ir/2, A = B 
and the LHS vanishes, so C — 0. Now set 



/ fa? a 2 + 1 /a; 2 a 2 + 1 

= csc~ 1 J — 1 _ a 2 — ' y=Vl-a~ 2 and u = J — - - ^ — . 
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da; 



o \/ x 2 + 1 — a 2 



: tan 



_j / V x 2 + 1 — a 2 



Ti 2 



V6 2 + 1 - a 2 + b 
1 + a 



-Ti 2 



V6 2 + 1 - a 2 - b 
1 + a 



(B.4) 



We next define 



g(a,b) 

Then g(a, 0) = (7r/2)Log(l + a) and 
5 



s 2 + 1 



Log[V& 2 (s 2 + l) + l + a]. 



96 



Sf(o,6) = b 



[ v / b 2 (s 2 + l) + l-a] du tan" 1 ^ 2 + I -a 2 /a) 



o (6 2 s 2 + b 2 + 1 - a 2 ) v /6 2 (a 2 + 1) + 1 V& 2 + 1 - a 2 



Reintegration using ( IB.4i gives 



7(o, ft) = -ln(l + o) + Ti a 



V6 2 + 1 - a 2 + b 
1 + a 



-Ti 2 



+ 1 - a 2 - b 
1 + a 



which is easily transformed into 



Logf-s/z 2 + c 2 + f i 



i da? = jLog[l + Vc 2 - a 2 ] 
.x +1 2 



Ti 2 



\/c 2 + a 



Vc 2 - a 2 + 1 



+ Ti 2 



1 + \/c 2 - a 2 
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